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Abstract. We characterize in details the aging properties of the ferroelectric phase of K Tai-a, Nb^ O3 
(ktn), where both rejuvenation and (partial) memory are observed. In particular, we carefully examine 
the frequency dependence of several quantities that characterize aging, rejuvenation and memory. We 
find a marked subaging behaviour, with an a.c. dielectric susceptiblity scaling as ujy/tZ, where i w is the 
waiting time. We suggest an interpretation in terms of pinned domain walls, much along the lines proposed 
for aging in a disordered ferromagnet, where both domain wall reconformations and overall (cumulative) 
domain growth are needed to rationalize the experimental findings. 

PACS. 05.40.-a Fluctuation phenomena, random processes, noise, and Brownian motion - 77-22-Gm Di- 
electric loss and relaxation - 78-30-Ly Disordered solids 



1 Introduction 

Aging is a widespread phenomenon! which manifests itself 
through the dependence of the properties of a material on 

\ the history of the studied sample [Q . It has been observed 
and widely documented on the magnetic susceptibility of 

' spin-glasses (SG) ||-|| and disordered ferromagnets Q, 
on the elastic compliance of polymers [Q , on the dielectric 
constant of disordered dielectrics ('dipolar glasses') J|,[5| 
and of some liquids p| , and recently on the rheological 
properties and dynamical structure factor of 'pasty' col- 
loids |I^Jl3| ]. Aging reflects the fact that the time needed 
for the system to equilibrate becomes larger than the ex- 
perimental time scales. Although several theoretical mod- 
els for aging have been investigated, the correct physical 
picture for aging is still rather controversial, in particular 
in spin-glasses. The simplest model for aging is that of do- 
main growth: below some critical temperature, the system 
progressively orders by growing domains of the competing 
low temperature phases (for example up and down in the 
case of an Ising ferromagnet). In this case, the aging of 
the susceptibility can be understood from the dynamics 
of pinned domain walls: as time elapses, the density of do- 
main walls decreases, and the remaining domain walls are 
better and better pinned. 

Stronger constraints on microscopic theories appear 
when one wants to interpret the so-called 'rejuvenation 
and memory' effect [Q. Suppose that one cools a glassy 



system at a constant cooling rate 1Z down to a certain tem- 
perature T p i, which is then held constant. At this tempera- 
ture, aging manifests itself by the decrease of some suscep- 
tibility (both of its real and imaginary part). If cooling is 
then resumed after a 'plateau time' i w , the susceptibility is 
seen to increase again, approaching the value it would have 
had if i w was zero: the system 'rejuvenates' as it seems to 
forget its stay at T p \. If cooling is carried further, down to 
low temperatures and then followed by a steady heating, 
a 'dip' in the susceptibility is observed when passing back 
through the temperature T p \. This shows that the system 
has actually kept some 'memory' of its evolution at T p \. 
The coexistence of these two effects is quite remarkable, 
and has led to several theoretical scenarios |l5Ul6| , [l7| , ^8| . 
These two effects have been seen in a variety of materials: 
various spin-glasses, pmma (plexiglass) (p]] or ktn (a dis- 
ordered dielectric, see below) [^0|. However, some materi- 
als exhibit aging but no rejuvenation or rejuvenation 
but no (or partial) memory 

The aim of this paper is to characterize in more de- 
tails the aging properties of the ferroelectric phase of KTN, 
where both rejuvenation and (partial) memory are ob- 
served |p0|| . In particular, we have carefully examined the 
frequency dependence of several quantities which charac- 
terize aging, rejuvenation and memory. We suggest an in- 
terpretation in terms of pinned domain walls, much along 
the lines proposed in || for a disordered ferromagnet, 
where both domain wall reconformations and overall do- 
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main growth are needed to rationalize the experimental 
findings. 



2 Experiments 



The pure potassium tantalate crystal K Ta O3 belongs 
to the cubic perovskite family. Quantum fluctuations pre- 
vents a ferroelectric order to establish fully, even at T = 
OK. If a fraction x of tantalum ions is randomly substi- 
tuted by isoelectronic niobium ions to give crystals of K 
Tai-z Nbz O3 (ktn), the trend towards ferroelectricity 
of the perovskite lattice is enhanced. It has been assumed 
that this tendency was due to random fields generated by 
the niobium ions [ p2[ . Indeed, if the niobium concentra- 
tion x is larger than x c — 0.008, the cooperative ordering 
gives way to a ferroelectric phase at temperatures below 
a finite transition temperature T c . For the sample used in 
the present study we observe a broad transition at T c = 31 
K, defined by the maxima of the real part e' and the imag- 
inary part e" of the dielectric constant, independently of 
the measuring angular frequency u>. When x > x c , there 
are actually several transitions: when the temperature is 
lowered, ktn is successively cubic, tetragonal, orthorhom- 
bic and finally, rhombohedral p3[ . For the concentration 
x of our sample, the three transitions merge into a single 
one. From the phase diagram obtained in j23|, we infer 
that the niobium concentration is around x = 0.022. The 
experiments reported below were all performed at tem- 
peratures T < 25 K. Therefore, they concern an ordered 
(ferroelectric) phase which contains some amount of dis- 
order (the randomly substituted niobium ions). 

Using a Hewlett-Packard 4192A impedance analyser, 
we have measured the electric capacitance and the dielec- 
tric loss at seven frequencies, ranging from / — 1 kHz to 
/ = 1 MHz. These can easily be transformed into the real 
part e' and the imaginary part e" of the complex dielectric 
constant. We measure the complex capacitance C(w, t, T) 
as a function of time t while the sample temperature T is 
a controled function of time. Several procedures were used 
after an annealing near 55 K and an initial rapid cooling 
across the transition temperature T c down to T max = 22 
K. The simplest procedure was composed of a cooling at 
a constant rate dT/dt = —TZ from T max to T m i n = 4.8K, 
immediately followed by a steady heating at the opposite 
rate dT/dt = +TZ from T m ; n to T max : this was used as the 
reference curves. Another procedure is to cool at the same 
rate —TZ between T max and an intermediate temperature 
T p i. There, cooling is interrupted by an isothermal evolu- 
tion or plateau at this temperature, which lasts typically 
t w = 10000 seconds. Cooling then resumes down to T min 
and is again immediately followed by a steady heating at 
the opposite rate from T m ; n to T max . A typical cooling rate 
is TZ — 5.2 mK/s. It was carefully checked that the linear 
response conditions are fulfilled. 
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Fig. 1. Behaviour of the total decay of the real-part of the 
capacitance, 8C'(u>, t w , Tpi) = C'(uj,t = + ,T p i) — C'(uj,t = 
iw,Tpi), as a function of the frequency / = uj/2h, for i w = 
20000 s and T p i = 12.1 K. The magnitude of this decay is 
found to behave as a power-law of the frequency: SC' oc ui~ v , 
with v = 0.127±0.005 over three decades in frequency (1 kHz - 
1 MHz.). We show in the inset the ratio SC" /SC 1 as a function 
of frequency. For a power-law dependence, this ratio should be 
constant, equal to tan(i/7r/2) = 0.202, shown as the full line. 



3 Frequency analysis of isothermal aging 

The isothermal variation of the complex capacitance C(w, t) 
in KTN is a slow decay in time (which also depends on the 
cooling rate TZ @). This decay can be fitted by a power 
law (to + t)~ v with a small exponent r\ (typically ~ 0.05) 
and to in the range of 10 2 s. For a given plateau temper- 
ature T p i and a given plateau duration t w , the total de- 
cay of the complex capacitance, 5C(u>, t w ,T p \) = C(u>, t = 
+ , Tpi) — C(u>, t — t w , T p i) > is found to obey a negative 
power law of the frequency / and therefore of u> = 2nf. 
The real part and the imaginary part of this variation can 
be fitted by (see Figure 1): 



8C'{u,U,T pX ) = N'uj- 



5C"(u;,t w ,T pl )=N"u;- 



with the same exponent v. This is actually a general con- 
sequence of the Kramers-Kronig relations, which hold for 
a linear, stationary and causal response theory. Although 
aging means that the response is non stationary, these 
Kramers-Kronig relations are valid in the quasi stationary 
limit where oj^ 1 <C t w , which obviously holds in our ex- 
periments (where tot w 3> 10 3 ). The Kramers-Kronig equa- 
tions actually also provide a precise relation between the 
amplitudes N' and N", which only depends on the value 
of v. 

^ = tan (t) ■ (2) 

This relation provides a second determination of v, by 
checking the constancy of the ratio of the imaginary and 
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Fig. 2. Rescaling of the capacitance for various frequencies u) 
(from 1 kHz to 1 MHz) and plateau times f w . The temperature 
here is T p \ = 9.8 K. We have chosen the stationary part C s t(u>) 
to obtain the best rescaling, as a function of (cu 1 ^' J 't w )~ v , with 
H = 1/2 and rj = 0.05, corresponding to v = 0.10. The rescaling 
is not perfect for the early time part where the time scale to 
introduced in the text cannot be neglected. 
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Fig. 3. Capacitance C'(ui,T) for two thermal histories, one 
without a plateau at T p \, the other with a long plateau at T p i. 
The important points A, B, F used in the text are defined 
in this figure. The inset shows a zoom of the region around T p \. 
Note the small temperature fluctuations during the plateau 
period. 



real parts of 5C, and using Eq. (g) to extract v. As shown 
in Fig. 1, obtained for T pl = 12.1 K and t w = 20000 s, 
these two procedures allow one to obtain rather close es- 
timates for v. The regression in log-scale based on Eq. (0) 
gives v — 0.127±0.005. The exponent v is found to depend 
weakly both on temperature and time. More precisely, it 
varies linearly with T p i but only logarithmically with t w . 
The magnitudes N' and N" are increasing functions of T p \ 
and t w . 

In spin glasses, the dependencies of the aging part of 
the susceptibility on the waiting time t w and on the fre- 
quency u> are not independent, but rather follow a simple 
ujt w scaling, which means that the typical relaxation time 
after waiting t w is t w itself. More generally, one may ex- 
pect scaling with u)t^, where \x is an exponent, generally 
found to be smaller than 1: this corresponds to 'sub-aging' 
|p|,P^. It is interesting to look for a similar scaling in our 
data. As shown in Fig. 2, the aging part of the capaci- 
tance can indeed be approximately rescaled with ^ ~ 0.5, 
leading to: 

C{uj,t w ,T pl ) ~ (w, Tpi) +C(T pl ) x (wt£)-" t w » t , 

where C s t is the stationary part of the response. This equa- 
tion describes simultaneously both the above discussed 
frequency dependence and the time relaxation of the ag- 
ing part, with the exponent relation rj — v[i. 



4 Frequency analysis of rejuvenation and 
memory 

A useful way to parametrize the 'strength' of rejuvenation 
is to compare the slope S — dC'/dT of the capacitance 
with temperature just after the plateau at T p i with the 
slope S r of the reference curve, which is positive for all 
frequencies. The stronger the difference SS = S r — S, the 
stronger the 'rejuvenation' effect. We find again that SS 
is a negative power-law of the frequency: small frequen- 
cies are more efficiently 'rejuvenated'. The value of the 
exponent is found to be similar to that of v. For example, 
for Tpi = 14.3 K and t w = 10000 s., we find an exponent 
equal to 0.20. Note that for small enough frequencies, the 
capacitance increases when the temperature is decreased 
(S < 0), while the capacitance decreases (S > 0) at high 
frequencies, as is the case of the reference curve. 

Figure 3 shows, as a function of temperature, the real 
part of the capacitance C'(ui, T) for two thermal histories. 
The first curve is the reference curve corresponding to a 
regular cooling from T max to T m - m , followed by a regular 
heating from T mm to T max . The value of C'(u>, T m - m ) along 
this path is indicated by point 'C in Figure 3, and the 
value of C'(u>, T pl ) at T = T p \ are point 'A' on the cooling 
curve and point '-E' on the heating part. 

The second curve corresponds to a regular cooling from 
Tmax to Tpi (point 'A' again), followed by isothermal ag- 
ing at T p i = 12.1 K for a certain time t w = 20000 s, 
reaching point '£?'. The isothermal decay reported above 
is therefore SC'(ui, £ w , T p i) = \A — B\. Then regular cool- 
ing resumes from T p i to T min , where it reaches point '£)' 
which is below point 'C. The distance |C — D\ reflects 
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Fig. 4. The memory indicator M'(u>) as a function of fre- 
quency for T p i = 12.1 K and t w = 20000s, in a log-log plot. 
The exponent found here is the same as in Fig. 1, as demon- 
strated in the inset where the 'memory ratio' M! /8C' is shown 
to be constant, around 0.5, over at least three decades in fre- 
quency. 



some cumulative aging due to the plateau at T p \. The last 
part of the curve is again obtained by a regular heating 
of the system. The value of C'(u), T p \) when passing T p \ is 
point 'F'. 

If the system had lost all specific memory of its stay 
at temperature T p i, point F would be situated at point E 
translated downwards by the cumulative part of the mem- 
ory, measured by | C — D \ . The extra difference can thus 
be taken as a measure of the specific memory of the pro- 
cesses taking place at T v \. We therefore define the memory 
indicator AA'{ui) as: 



M'{u) = \E-F\- \C-D\. 



(4) 



A situation where the memory is perfectly conserved would 
lead to M 1 = \A- B\. 

As shown in Fig. 4, this difference is positive in KTN, 
and again behaves as a power-law with frequency. Interest- 
ingly, the exponent describing this behaviour is very close 
to v. For T p i = 12.1 K and t w = 20000 s, we find a value 
of 0.13, nearly identical to the value of v reported in Fig. 
1. In the inset of Fig. 4, we have shown the 'memory ra- 
tio', i.e. the ratio of M! to the total isothermal decay 5C 
during the plateau studied above, for different frequen- 
cies. This ratio is around 0.5 (showing that some memory 
is lost), but is nearly independent of frequencies, showing 
that the frequency dependence of the memory is actually 
that of the aging part of C . A simple interpretation of 
this fact will be given below. 



5 Interpretation 

The microscopic origin of aging in disordered ferromag- 
nets or ferroelectrics is most probably the slow dynamics 
of domain walls, pinned by impurities. The important in- 
gredient to understand the physics of these pinned objects 
is the fact that the typical pinning energy associated with 
a portion of wall of linear length £ grows as Y£® where 
9 is a certain exponent and T an energy scale Ea|. This 
means that the time needed to equilibrate such a portion 
of wall grows exponentially, as t(£) ~ exp(Y£ e /T). Corre- 
spondingly, as emphasized in |pl| , time scales are strongly 
hierarchized, in the sense that t(2£) S> t(£). Hence, some 
small scale reconformations can take place while the over- 
all, large scale configuration of the domain wall is fixed. 
This small scale motion between metastable states leads 
to better and better pinned configurations, and therefore 
a smaller susceptibility: this leads to aging, even if the 
overall size of the domains does not change. The large 
scale coarsening leads to a progressive increase of the av- 
erage domain size R, and thus a decrease of the fraction 
of dipoles belonging to domain walls, as 1/R. The coexis- 
tence of these two effects suggests that the aging part of 
the susceptibility (magnetic or dielectric) should read: 



Cag(w, t w , Tpi) 



1 



t(w,Tpi)+c rec (u;,t w ,Tpi)], (5) 



where c st (w,T p i) is the equilibrium contribution of a wall 
of fixed (large) size, and c iec (u),t w ,T p i) is the aging con- 
tribution, due to small scale reconformations which al- 
low the wall to find deeper and deeper metastable states. 
Simple 'trap' models for these reconformations suggest 
that c rcc oc (u>t w )~ 1 ' {p < l)|p6|. More general (subag- 
ing) forms are possible fl(|^7j , as observed experimentally 
(see Figure 2). A possible physical mechanism is that each 
metastable state is visited a large number of times. Note 
that in Fig. 2 we have neglected the time dependence of R, 
which is justified if the cooling rate \TZ\ is small enough. 
In that case, the relative growth of R during £ w can be 
neglected. 

As explained in 26 16| , |l8| ], the reconformation con- 
tribution can account for the rejuvenation and memory 
effect. As the temperature is decreased, the small scale 
motion, which was in equilibrium at the higher tempera- 
ture, is suddenly driven out of equilibrium. The aging part 
c roc (ujt w )~ u has a singularity for t w — that reflects 
the excess dissipation that arises in highly out of equilib- 
rium situations. Therefore any 'micro-quench' leads to an 
increased susceptibility and new aging dynamics. Simul- 
taneously, larger scale dynamics is effectively frozen and 
contributes to memory. The distinctive feature of Eq. (|^) 
is however the cumulative factor 1/R that always grows as 
soon as the system is in its ferroelectric phase (albeit more 
slowly at smaller temperatures) . This cumulative factor is 
at the origin of all cooling rate effects (see the discussion 
in |2]| ) . It also explains the change of sign of the slope S 
reported above: at high frequency, the reconformation con- 
tribution c rcc (co>, t w , Tpi) is negligible, and the main effect 
is the progressive increase of R which leads to a decrease 
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of C ag — c st (iLJ,Tpi)/R(t). Therefore the slope S just after 
the plateau is positive: the capacitance decreases when the 
temperature is lowered further. On the contrary, the reju- 
venation effect is dominant at low enough frequencies, and 
explains why S is negative. The frequency dependence of 
5S is thus expected to be similar to that of the isothermal 
decay of the capacitance, as was reported in section 4. 

It is important to note that the large scale domain mo- 
tion is also responsible for the loss of memory, as suggested 
in ||: if, during the period spent at lower temperatures, 
the domains have had enough time to move substantially 
and lose track of their previous position relative to the 
impurities, the memory effect is lost More precisely, let us 
call p(t\ ow ) the probability that a given wall has not sub- 
stantially moved during the time ii ow = 2(T p i — T m i n )/|7?.| 
spent at low temperatures, and thus retains the memory 
of the past history. The capacitance reached at T v \ on the 
heating curve after a plateau at T p \ (point 'F' in Fig. 3) 
is therefore equal to: 



Cpj(w) 



1 



R{tw + £iow) 
+ (1 -p)c iec (uj,Q + ,T pl 



(w, t w , T p i) + 



(0) 



where the origin of time is chosen when T p \ is first reached 
(point 'A'). The above equation means that a fraction p 
of the signal has kept full memory, while a fraction 1 — p is 
totally rejuvenated. This decomposition is only approxi- 
mate, since the memory part is affected by the time spent 
near T p \. The memory indicator defined above can be writ- 
ten as: 



M(u) = 



1 



R(t w + iiow) 



[c st (w,T pl ) 



+c roc (u>,0 + ,Xpi)] - C p i(u), 



(7) 



where the first term in the right hand side is the capaci- 
tance that one would measure in the absence of any mem- 
ory (p = 0). We there find that the memory indicator 
M is related to the isothermal decay of the capacitance 
6C(uj,t w ,T p i) by: 



M{w) 



p(iiow) 



R{t w + tl ow ) 



[c roc (w, + , T p i) - c rcc (w, t w , T pl )] , 

(8) 



or 



M(u) = p(tlov,) 



R(t v 



R(t w + <low) 



5C(cj,t w ,T pl ). (9) 



This proportionality is indeed confirmed by Fig. 4. Note 
that p is expected to decrease when ti ow increases. On the 
other hand, the ratio i?(i w )/i?(i w -I- t\ ow ) (which is ex- 
pected to be close to 1 in te experiments) also decreases 
when t\ ow increases. Therefore, the amplitude of the mem- 
ory is expected to decrease as the cooling rate decreases. 
We have checked this feature experimentally, and have 
found that the dependence of on the cooling rate 1Z 
can be fitted approximately by a power law M.' cx TZr with 
(3 = 0.28 for T p i = 14.4 K and t w = 10000 s., and TZ in the 
range 1.3 to 20.8 mK/s. This experimental dependence 



suggests that for very small cooling rates, the memory 
effect is eventually totally erased by domain growth, as 
proposed in [fj for disordered ferromagnets. 



6 Conclusion 

The present experiments can be seen as a quantitative 
confirmation of the qualitative scenario of aging proposed 
in [|||l6| for disordered ferromagnets, which ascribes the 
observed effects to the motion of pinned domain walls. 
The simultaneous presence of overall domain growth and 
internal reconformation modes allows one to account for 
the phenomenology of temperature cycling: cumulative ef- 
fects and memory erosion are due to the former, whereas 
rejuvenation and memory are due to the latter, and there- 
fore show up for small frequencies. We have also shown 
that all the quantities that measure the reconformation 
contribution behave as power-laws of frequency, and that 
a subaging (wi(£) scaling holds, with a rather small value 
of the exponent p ~ 0.5. 
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